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Phase-Sensitive Impurity Effects in Vortex Core 
of Moderately Clean Chiral Superconductors 
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We study impurity effects in vortex core of two-dimensional moderately clean superconductors 
within the quasiclassical theory. The impurity scattering rate T(E) of the Andreev bound states 
in vortex core with +1 vorticity of p-wave superconductors with d = z (p x + ip y ) is suppressed, 
compared to the normal state scattering rate F n in the energy region F^/Eg <C E <C Eg = |5o|Aoo 
with scattering phase shift So (\5o\ -C 1) and the pair-potential in bulk Aoo. Further we find 
that T(E)/T n for p-wave superconductors with d — z (p x — \p y ) is at most 0(£'/A 00 ). These 
results are in marked contrast to the even-parity case (s,d-wave), where T(E)/T n is known 
to be proportional to ln(Aao/E) . Parity- and chirality-dependences of impurity effects are 
attributed to the Andreev reflections involved in the impurity-induced scattering between bound 
states. Implications for the flux flow conductivity is also discussed. Novel enhancement of flux 
flow conductivity is expected to occur at T <C E$ for d — z (p x +ip y ) and at T <C A^, for 
d - z (p x - ip y ). 

KEYWORDS: Superconductivity, vortex, impurity effect, quasiclassical theory, flux flow conductivity 
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§1. Introduction 

Low energy physics in vortex state of superconductors 
without gap nodes is governed by localized excitations 
inside vortex cores. The density of states within a core 
is comparable to that of the normal phasalP and hence 
vortex cores may be regarded as locally realized normal 
region. This picture has been used in the calculation of 
the flux flow conductivity in ref. 0. Strictly speaking, 
however, the description of vortex core as a normal re- 
gion is limited to the dirty superconductors, where the 
coherence length £o is much larger than the mean free 
path I. 

In clean superconductors where £,o < I is satisfied, par- 
ticles and holes within a vortex core are subject to An- 
dreev reflections and the coherent superposition of par- 
ticle and hole states form bound states before they are 
scattered by an impurity. The flux flow conductivity erf 
for clean s-wave superconductors has been calculated in 
refs. |[[|,|] and found to be 



<jf /ct„ ~ In (Aoo/T) H c2 /B, 



(1.1) 



in terms of the normal conductivity a n , the upper criti- 
cal field H C 2 and magnetic induction B. The logarithmic 
factor in rtO) results from the shrinkage of vortex core at 
low temperatures (Kramer-Pesch effect)© and logarith- 
mic energy dependence of the impurity scattering rate of 
the Andreev bound states.utP 

Several clean superconductors are realized in corre- 
lated electron systems, such as UPt3, organic conduc- 
tors, cuprates and Sr2Ru0 4 . Those systems are_believed 
to have unconventional superconducting statesJj' 

Flux flow conductivity of d-wave superconductors has 
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been calculated in refs. [T^,|T^ and it is similar to the 
s-wave result (1.1) in moderately clean regime, where 



ZAoo/ep < £ < I, 



(ep; Fermi energy) 



(1.2) 



is satisfied. Thus the flux flow conductivity of even- 
parity superconductors are expected to be similar. 

However, flux flow conductivity of odd-parity super- 
conductors may be different. Actually, it has been shown 
in ref. [L2 that single impurity does not change the spec- 
trum within a vortex of a chiral superconductor when 
only the backward scattering is considered. Further, 
from the random-matrix approach, level statistics within 
a p-wave vortex core is expected to be different from 
that within an s-wave coreO'Lj Even experimentally, 



appreciable deviation from (1.1) has been reported in 
the measurement in UPt 3 , where odd -pa rity supercon- 
ducting states are presumably realized.E£IEj' Thus novel 
impurity effects in vortex core and flux flow conductivity 
are expected in odd-parity superconductors 

For lastiive years, Kopnin et a/.,Lrl3' StoneEZP and Blat- 
ter et alx^ developed the kinetic theory for vortex dy- 
namics in clean superconductors. Their results on flux 
flow conductivity in the moderately clean regime reduce 
to 



<7f ~ neojQT/B, 



(1.3) 



with electron density n, charge unit e(> 0) and the re- 
laxation time t and angular velocity lu$ of precession of 
the Andreev bound states within vortex cores x3 Since 
uj is expected to be of the order of In (A m /T) Aj^/ep, 
the calculation of tJf, essentially, reduces to that of re- 
laxation time t; impurity nmblem. inside vortex cores. 

Although several authorsEroliiJ addressed the impu- 
rity problem within vortex of odd-parity superconduc- 
tors, microscopic theory valid in the moderately clean 
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regime (1.2) is still rrj.issing.c3 1 Such a theory is highly 
desirable in the following reasons: fhjst, an existing odd- 
parity superconductor Sr2Ru04 hasc2P Aoo/ep ~ 10~ 4 
and Z/^o ~ 8 and thus belongs to the moderately clean 
regime. Another odd-parity one UPt3 hascJ A^/ep ~ 
0.05 and l/^o ~ 50 and lies between (2) and superclean 
regime (I 3> £o e F/Aoo). Study in the moderately clean 
regime, therefore, has relevance to experiments. Sec- 
ond, it is not so clear to what extent the approxima- 
tion adopted in ref. Ill is valid—. Third, the validity of 
the random-matrix approacriiS E3' can be examined only 
from microscopic theory. 

In the present paper, we consider the impurity prob- 
lem within pancake vortex cores of two-dimensional 
s-wave and chiral p-wave superconductors with d = 
z (p x ± ipffl wit, hi n the quasiclassical theory of supercon- 
ductivityOa'Oo' Although the method to calculate 
the impurity scattering rate has been given in refs. p}|26[ 
we develop an alternative analytical method, which is 
based, pa—B.icfiati— formulation of the quasiclassical the- 
ory B-ByHy bk3> 

In the next section, we present the framework of qua- 
siclassical theory of superconductivity. In §3, impurity 
problem of s-wave vortex core is considered. In §4, the 
results on p-wave case are given. In §5, a brief argu- 
ment on the flux flow conductivity is given. In §6, we 
discuss physical origin of the parity- and chirality- de- 
pendent impurity scattering rate inside vortex core. We 
also discuss the formulation used in §3 and §4. In §7, we 
summarize the present study. 

§2. Quasiclassical Theory of Superconductivity 

We consider two-dimensional superconductors with 
circular Fermi surface in Uma_U_ limit. Quasiclassical 
theory of superconductivitycj'SoL^P is described in the 
equilibrium case by the quasiclassical green function 



9(z,P,r) 



f) 
■f 



f 
-9 



(2.1) 



which is a 2 x 2 matrix in particle-hole space and is a func- 
tion of (complex) frequency z, direction p = (cos a, sin a) 
of momentum p — pfP, a point r = r(cos (f>, sin <j>) in real 
space. The equation of motion for g is given by 



- iv(p) ■ Vg = |er 3 - A(r,p) - £,g| , 
supplemented by the normalization condition 

.9 = -71- 1, 

where 1 is the 2x2 unit matrix. 
Here T3 is a Pauli matrix, 

1 
-1 



and A is given by 



A(r,p) 
-A*(r,p) 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



where A(r,p) is the pair-potential. The impurity self- 



energy S in the t-matrix approximation is given by 



£(z,r) 



1 



(2.6) 



N V(g) ' 

with impurity concentration m, the density of states in 
the normal phase No and the scattering potential for a 
single electron V. The symbol (• • •) denotes the average 
of ■ • • over the Fermi surface. 

For later convenience, we parameterize the self-energy 



(2.6) in terms of impurity scattering rate of normal phase 

1 riiN a V 2 



r n = 



2ttt„ 



1 



(2.7) 



(TrNoVy 

and t he s cattering phase shift tan5 = —ttN V. Expres- 
sion (2J3) then becomes 

6 r n(.9) 



COS 2 £q 



'■So ((f)(1) -(g)'- 



(2.8) 



where we have dropped the term proportional to 1; it 
gives no contribution in (2.2). In the Born limit (<5 — ► 0), 



the self-energy reduces to 

t = T u (g). 



(2.9) 



In this paper, we consider the case where g is an ana- 
lytic function of z in the upper half complex plane. Set- 
ting z = e + i<5 with real e, we obtain the retarded Green 
function. 

In our case, the Fermi surface is circular and hence, 
the Fermi velocity v(p) is given by vp with a constant v. 
The left hand side of ( |2.2] ) then becomes —wdg/ds with 



s = r ■ p = r cos (<^> — a) . 



(2.10) 



Thus, the direction of_momentum specifies a line ("qua- 
siclassical trajectory"! 2 ^), on which the quasiclassical 
equation is to be solved. Each quasiclassical trajectory 
can be determined by the condition b(r,a) — constant, 
where the impact parameter b(r, a) is defined as 

b(r, a) — r ■ (z x p) = r sin (<j> — a) . (2-11) 

In this paper, we solve the quasiclassical equation 

ion of the quasiclas- 
'■' The solution g of 



(2.2) using a special, 
sical Green function.! 



(2.2) can be written as 

1-77 
2 7 t 



1 



-77 



27 



-1 



77' 



(2.12) 



Here 7 and 7T are the solutions of the following Riccati 
differential equations: 

iv ■ V 7 = - (A* - E 21 ) 7 2 - 2 (e - En) 7 - (A + £y)3a) 

iv • V 7 f = - (A + E 12 ) 7 t2 + 2 (e - E n ) T * - (A* ^2513J>) 

In this parameterization, the normalization condition 



(2.3) is automatically satisfied. Now we fix e, a and 
b = 7'sin( — a ). Along the quasiclassical trajectory, 
equations ( 2.13| a,b) are to be solved respectively, under 
the initial conditions: 
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vt = 



-e-i5-i\/\A\-e 2 
A* 



for s = -(2 c 14a) 
for s = s c ,(2.14b) 



e + i<S + i\/|A|-e 



where s c is a positive number much larger than £o and 
eventually t he li mit s c — > +00 will be taken. The initial 
co nditi ons ( 2.14| a,b) are obtained from the bulk solution 
of (2.2). In those initial conditions, the effect of nonmag- 
netic impurity is absent as a result of Abrikosov-Gorkov 
theory for isotropic s-wave homogeneous superconduc- 
tor s.O 

We consider a single vortex with +1 vorticity using 
the following test pair-potentials: 



A(r,p) 

AoMe 1 ' 



1 for s-wave superconductors(SC). _^ 

e ±i« for p-wave SC (d = z(p x ± ifiy)) > 



The amplitude of the pair-potential Ao(r) is assumed to 
be a function of r = |r| (axisymmctric) and behave as 



A (r) = 



Aoor/^i, forr<Ci 
Aoo, for r > £1 



(2.16) 



with a characteristic length £1 and the amplitude of the 
pair-potential in the bulk A^. We take £i(< £0) as a 
length different from the conventional coherence length 
£0 = w/^Aoo) in general, in order to take account of the 
shrinkage of a vortex corefl^ 

§3. Isotropic S-wave Case 

3.1 Self-energy in Born limit 

In pure superconductors (S = Ol,|-thc approximate 
expression for g has been obtained.BE-J' For e <C A m 
and b <C £0 , 

5o(e,a,s,6) 



.9 



e + i<J - E(b) 



+ (regular part), 



(3.1) 



Here the spectral function go{e, a, s, b) has no singularity 
and varies with energy over the range of Aoo and with 
s over that of £o- E(b) is the pole corresponding to the 
Andreev bound state for a given quasiclassical trajectory. 



E P (b) 



with 



and the normalization 



c 


/o s 


2feA ~ m 


V2£i 




2 Z 1 ! 8 ! 
u(s) = — / ds'A( 
u Jo 


>(*') 





c 



dsexp[— u(s)] ~ tt£o/2 



(3.2) 



(3.3) 



(3.4) 



(the subscript "p" in (p.2| ) stands for pure superconduc- 
tors). The second term in the right hand side of (3.1) 
comes from the scattering state of quasiparticle. This 
contribution can be neglected in the case of e <C A^. 

Next we consider the case with impurity. For 
Max (e, r) <C A^, main contribution to g comes from the 



Andreev bound states and an expression similar to (3.1) 
still holds. In this case, however, those bound states are 
subject to the impurity scattering and the bound state 
pole moves to lower complex plane. Now we assume that 
g for Max (e, T n ) <C A M and b <^C £0 has a form of 

A g {e,a,s,b;T n ) 

9 ~ e + iT(b)-E(b)' [6 - b) 

with the weight function go varies with e or T n over the 
range of A^ and with b over the range of £q- The im- 
purity scattering rate T(b) for Andreev bound states ap- 
pears in (3.5) as the imaginary part of pole of Green 
function. 

Now we calculate the self-energy S a s a function of 



(e, a, r; T n ) by taking the average of ( 3.5) wit h resp ect t o 



a for fixed r. s, b depend on a through ( 2.10 ) and ( 2.11 ), 
respectively. The weight function go can be replaced 
by the value for e = b = T n = with the accuracy 
of (9(e/A 00 ,r n /A 00 ,&/£o)- ThJS-.value is available from 
earlier results on pure case.ErlZrE}' The self-energy can 
then be written as 



S(r;T n ) 



7rtT n exp [— u(s)] Mo 
2C \e-E(b) + iT(b)) 



>, 



with ( pT0| ), ( pill ) and 
M = Mo (a) = 



1 



(3.6) 



(3.7) 



From now on, we assume that 
E(b) ~E + E'b, E'v/A^ = 



with 6-independent Eo and E' and 

r(6) < e 1 . 



A 00 >Max(|e|,r). 

(3.8) 



(3.9) 

The denominator in ( |3.6D leads to a distribution local- 
ized around a — a\, CX2 satisfying 



E(b = rsm((f>-a l )) i = 1,2. 



(3.10) 



Therefore "the slowly varying part "exp[— u(s)] can be 
replaced by 

exp [-u(si,2 = rcos(4> - 01,2))] ■ 



If we set bi = r sin(0 — a), then, (3.10) leads to b\ = b? 
and si = —S2- It follows that u(si) = u(s2) because u(s) 
i s an even function of s by definition ( |3.3| ) . The solutions 
( 3.10 ), however, exist only if r > (Ree — E ) /E'. As a 



result, the self-energy (3J3) becomes 

7rur n exp [-u(si)] , M (a) 



Z(e,r;T n )r 
with e = e — Eo 



2C 



e- E'b' 



(3.11) 



ir and 
§1 = Max [si, (Ree - E ) / E 1 ] 



A direct calculation shows that 

(l/(£-E'b)) = 



E'ry/1 - w 
with w = e/(E'r). Here the branch of 



(3.12) 
is taken such 



1 
that 
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Similarly, we obtain 



—iw, 
1, 



for \w\ > 1 
for \w\ <C 1 



(3.13) 



kept. The quantities (3.1£) are given by 13 " 

d £ 7 ( .« = e"( s ) 






.±ia 



/(e-E'b + iT)) = 



=Fie 



±W 



^'rVl 



vr 



and 

9 h 7 ( ' f) 



3 u(«) 



(3.14) 
with the branch (3.13). 

By substituting (|3.11|) with ( p^ ) and Q3.14Q into J2^ ) 
and solving the resulting equation, we can obtain Eq, E' 
andT. 



-srV"^ 



T2V- 1 



T-s<: 



(3.20) 

ds's'- 1 A(s')e- u(s ' ) 
(3.21) 



In ( 3.2C ) and ( 3.21 ), the upper signs are respectively for 
<9 e 7 or <9b7 and the lower one for d e j' or dbrf ■ 



,-,( .t) 



Now we consider 7j» . They are the solutions of 

•2E, (3.22) 



as 



2E, V ^H = 2A7+ 



<9.s 



3.1? Bound state pole in Born limit 

Now we solve quasiclassical equation via (ETjra,b) with with 
a particular interest of the bound state pole. From 
(2.13a,b), the equation 

d /■„ +\ r, . * „x ,».„■> +-1 /- . -^Equations (3.22) are supplemented by the initial condi- 

tions 



E = E u - i (E12 + E21) /2. 



(3.23) 



-i^|-(l + 77 t ) = {(A*-S 21 )7+(A + E 12 )7t}(l+77tJ c 
us tic 

follows. This equation shows that the denominator of 
g vanishes for all s on the trajectory if it vanishes at 



7s(s = -s c ) = 0, 7|,(s = s c ) = 0. 



(3.24) 



a point S3 Thus the energy e which gives (I + 77TJ = Equations ( 3.22 ) with ( 3.24 ) yield the following solu- 
has the meaning of a bound state pole even in the tions: 
presence of impurity. 



What we need are the solutions for ( 2.13 a,b) in case 
of e <C Aoo, b <C £0 and E <C Aoo. The last inequality is 
equivalent to 



7^ 



= F 2w- 1 e tl(s) f ds'e- u( > s '^(s r ) 



(3.25) 



T„/ |e| « 1 



(3.15) 



because 



7TTn 

2C\e\ 



After t aking the l imit s c 
(|3~20| ), ( ^2l|) and ( ^25|) 

1 + 77 T = 1 + 77 
4Ce"( s ) 



00, we can obtain from 



C~£o 



7rA c 



1 



We expand equations ( 2.13 a,b) with e, b and E up to 
first order. It is convenient to introduce the following 
notations: 

7 = 7 e iQ , 7t= 7 t e - iQ! , A = Ae ia , 

A* = At e - iQ , Eu = S 12 e iQ , E 21 = E 21 e" iQ . 

(3.16) obtain 
First we consider the case where e = 0, 6 = and E = 
0. In this case, the solutions for ( 2.13j a,b) with initial 
conditions j(s — —s c ) = — i and ^(s = s c ) = — i are 
given by 



Ep{b) ~~ijJ ds ^( s K u(s) ) 3 - 26 ) 

In ( 3.26J ), the second term in the curly bracket has been 
given by (pi). 



Now we calculate the third term in {• • •} in ( 3.26 ). 
First we note that the off-diagonal elements Ei 2 + E21 
vanishes on the trajectory with b = and E ca n be 
replaced by En. With the use of ( 3.11 ) and ( 3. 12] ), we 



7 = 7 = — i. 



1 

7r^r n 

2iC 2 



d s E 11 (s)e- u(s) 

p -2u(s) 

ds- 



^/Ws* - e' 



(3.27) 



_ , , ... , . , , e^ — c^. 

(3.17) 

For s ~ £o, u(s) ~ 2/S.aos/v. Hence, the above expression 
Next we consider the first order corrections with respect can b e reduced to 
to e, b and E. 7 and 7^ are expanded as 



7 ( ' t) = -i + e (a £ 7 ( ' t) )+^(47 ( ' t) )+7^ t) 

+ C(e 2 ,6 2 ,E 2 ,e6,eE,6E), (3.18) 



2iC 2 J 



where 



a £7 ( ' f) = 



d 7 ( .t) 



dc 



a 67 ( A) = 



9 7 ( -t) 



ilh 



= r n 

with r n 



i 2e 



d s (£' 2 s 2 -e 2 ) _1/2 

_j Ime 



cot 



.0; ? 

Ree V A 2 



(3.28) 



7rwr n /(2C 2 J B')- From (iH)> (H3) and 
(3.28), we immediately see that 1 + 77^ vanishes at 



Ree = E p + nf n /2, Ime = -f n In 



2E V 



-r. (3.29) 



=b=S=0 "" e= fc=S=0 

(3.19) 
Here 7^ >" is either 7 or 7^. The symbol 7^ denotes 
the solution linear order of E in the following situation: 
b is taken to be zero, e as a coefficient of j( -t) is set Now we examine the consistency of our results. The 
to be zero while the implicit energy dependence of E is 



assumption ( |3.8[ ) is obviously satisfied in ( |3.29 ). 

The quantity T n is of the order of r n /m(7r£o/(2£i)) 
Therefore, the condition T <C A^ gives 
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M+ = M+ (a) = 



T/Aoo ~ (r n /AooJ lntAoo/Ep) 

~ (r n /A oo )ln(6/e )/ln( 7 re /(2ei)) « 53.30) 



Our assumption (3.15) on the self-energy can be rewrit- 
ten as 



r n /Max(£ p ,|e|)«l. 
The condition ( |3.9| ) is thus confirmed by 



(3.31) 



i_dr 

E''db 



E'b 



<1. 



£ p ln(7r£o/(2£i)) 

Now we can write down the final expression for g. 
The matrix part of g in (2.12) has a non- vanishing 
value for e = b = £ = 



1 — 77T 
2 7 t 



2l 



vt 



-1 + 77 
As a result, we obtain 

7rwe-" (s) M (a) 



2M (a) 



.9 



2C (e - E p - 7rf n /2 + IT 



(3.32) 



(3.33) 



under the condi tions (e,r n ) <C Aoo, fe <§; £0 and two 
additional ones ( 3.30J ) and ( |3.31 ). 



§4. P-wave Case 

In the present section, we consider p-wave supercon- 
ductors with the pair potential 

A(r) = A (r)e iW±Q) 

in a way similar to the previous case. 



(4.1) 



4-1 {Px + tyy) superconductors without impurities 

First we consider the green function in the pure case 
(£ = 0). The spectrum within a core of pure super- 
conductors with d = z(p x ± \py) has been calculated 
in ref. Bm numerically using the Bogoliubov-de Gennes 
equation. The result in the present subsection 4-1 corre- 
sponds to the /cfCo 3> 1 limit of that in ref. |35| . 

If we use the same parameterization (2.12) for g and 
the expression (4.1) and set £ = 0, equations ( 2.13| a,b) 
still hold. Further when we introduce the following no- 
tations: 



7 = 7c 



A* 



At e - i2Q , 



7 t = 7t c -i2« ; 



A = Ae i: 



(4.2) 



the resulting 7^ <*> in (4.2), then, turn out be exactly 



same as those in the previous section. We can thus use 



(3.20) and (3.21) to obtain the solutions for small e and 



small b. As a result, the expression for the quasiclassical 
Green of pure p-wave superconductors is given by 

nve- u ^M+ 



f) 



2C (e-E p + 16) ' 



(4.3) 



-i2a 



j2a 



1. 



(4.4) 



The fact that the phase factor in the off-diagonal part in 
g is not ±a but ±2a leads to a new aspect in the present 
case. 

4-2 (p x + ip y ) superconductors in the Born limit 

Now we consider the case with the self-energy We 
assume again that g has approximate form of 

9 ~ 2C (e - E Q - E'b + ir) ' ( ' } 

Following the same procedure as the previous section, we 
obtain from (4.5) the expressions for self-energy 

e- i2 *£ 12 = e i2 *£ 21 
Tive- U ^ ( \~2w 2 



2CE'r \^T^P 



2\w 



(4.6) 



For the diagonal part (En, £22), the expressions are the 
same as those for s-wave. 

Now we solve the quasiclassical equation with use of 
parameterization (2.12) ; 

-\tt ( 1 - 77 1 27e i2Q 
1 + W 7 V 2j^~ i2a -1 + 77 1 



9 = 



(4.7) 



7 and 7^ in (4.7) are the solutions of the following Riccati 
differential equations: 



iv ■ V7 
it; • V7' 



(At - E21) f - 2 (e - En) 7 - (A + E{4).8a) 
(A + S 12 ) 7^ + 2 (e - En) 7 f - (A f - pbg» 



where A and A* have been defined in (4.5). E12 and £21 
are defined as 



£i 2 e 



i2n 



E21 — E2ie 



-\2a 



(4.9) 



Equations (4.8a,b) should be solved under the initial con- 
ditions: 



7(s = -s c ) 



7 f (s = s c ) 



A(a = - Sc )/ c + iV|A|--{*W 



At( s = Sc )/ e - + iJ|A| 2 -e 2 (Ll0b) 



where e is given as the solution of 



e l-TriyVIAf-e 



(4.11) 



The initial conditions ( 4.10] a,b) are derived from the bulk 
solution for p-wave superconductors. Nonmagnetic im- 
purities give pair-breaking effects on the non-swave ho- 
mogeneous superconductors. The presence of impurities, 
therefore, affects the bulk solution through e. 

What we need are the solutions of (|4.8| a,b) for small e, 
b and T n . The initial conditions (4.10|a,b) are expanded 



7(s 



-s c ) ~ 7 f ( s = s c ) 
e 



o (er n , 



,(4.12) 
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The expression (4.12) shows that initial conditions for 
impure superconductors are the same as those for pure 
superconductors in the case of small e, b and T n ; this 
fact suggests that impurity effects should be different 
from those for bulk superconductors. 

Now we calculate 7 and y< for p-wave superconductors. 
All the results from ( p. 18 ) to ( 3.26| ) apply to the present 
case and hence we can begin with the expression ( 3.26J ) 
for 1 + 77T . 

From ( |4.6| ) for £12 and £21, it follows that £12 = £21 
on the trajectory with 6 = 0. As a result, those off- 
diagonal elements contribute to the bound state pole in 
p-wave case, in contrast to the s-wave case. 

The integral 



-1 
~C 



ds£(s)e 



-u(s) 



for p-wave case can be reduced to 



C 



d.s 



^£V2 S 2 _ g2 ^ E >2 S 2 _ g 2 _ ig ) 



(4.13) 



(4.14) 



The denominator of the integrand in (4.14) makes the 
integrand rapidly convergent. Therefore we can replace 



e 2u(s) ky un j^y After performing integration, (4.14) 
becomes i2r n . Consequently, we obtain 



E = 0, E' = Bp/b, r = 2f„ 



(4.15) 



for the green func tion (|4.5|) of p-wa ve superconductors. 
If we compare ( 4.15 ) with ( 3.29J ) , we can easily see a 
qualitative difference of T between the two cases; T in s- 
wave case has logarithmic dependence on the bound state 
energy while that in p-wave case is energy- independent. 
This difference comes from different behaviour of £ in 
(|3.26| ) and J4.14J) . In s-wave case, £(s) = £ s (s) ~ 1/s for 
|e| /E' <C s -c Co- This long tail gives the logarithmic de- 
pendence ln(£o-E"/ |e|) in I\ In £(s) of p-wave supercon- 
ductors, on the other hand, the long-tail of £n(s) is can- 
celled by that of £12. This leads to energy- independent 

r. 



4-3 (p x +ipy) superconductors away from the Born limit 
In this subsection, we consider the imp urit y scattering 
rate within the t-matrix approximation (2.8); the phase 
shift So is to be small (\6q\ <§; 1) but finite. Following 
the same method as that in the prev ious one, we obtain 
the green function in the form of (4J3) with Eq = 0, 
E' = E p /b and 



r 



(2ir/3 3 / 2 ) (E P /E S ) 



2/3 



for E s « E p < Aoo 
for T 3 JE$ « E p « E s , 

(4.16) 
We will describe the 



with E s = v\6 \/C ~ \5 \A X . 
derivation of ( j4.16| ) in Appendix. 

We see that finite phase shift sets an energy scale Eg, 
above which the Born approximation works well and be- 
low which it breaks down. The suppression of the impu- 
rity scattering rate in the energy region E p <C Eg may be 
related to "the spectrum rigidity" within a vortex core 
chiral superconductors against a single impurityjIiiJ The 
approximation in ref. n2l where only the backward scat- 
tering is taken into account, seems to be valid in that 



energy region. 

On the other hand, in the energy region Eg <C 
E p <C Aoo, the impurity scattering is large; this sug- 
gests the strong level mixing-arid, hence we expect that 
the random-matrix approachtJ'E3 should be applicable 
in this energy region. 

4-4 (Px — ipy) superconductors 

We consider in this subsection the impurity effects in 
(p x — ipy) superconductors under the test pair-potential 

A(r) = A (r)e^- a ». (4.17) 

We will see that the impurity effect in vortex core of su- 
perconductors with time-reversal-symmetry-breaking is 
chirality-sensitive. 

For those superconductors without impurity, we ob- 
tain in the same method for the green function 



nve 



- U ( S )M_ 



f] 



2C (e - E p + id) ' 



with 



M_ = 



1 -i 
-i -1. 



(4.18) 



(4.19) 



We note that the off-diagonal elements of the matrix part 
in (4.18) have no a-dependence. This fact will turn out 
to be crucial. 



On the basis of ( 4.18 ), we expect the green function of 
the form 

-< S )M_ 



9 ~ 2C(e-Eb-^ + ir) (4 - 20) 

in the presence of impurities. From ( 4.20| ), we immedi- 
ately see that 



Z-i\2 — ^21 — — iE 



11- 



(4.21) 



By following the argument from (3.18) to ( 3.26| ) with 
replacement of (7^ '" , £) by (7' '^, £) , we obtain for Eq 
and T in ( p~2C| ) 



£o=0, r = -/ ds£(s)e-"«. 



(4.22) 



The relation (4.21), however, yields £ = and accord- 
ingly, r in (4.22) vanishes. Our calculation is based on 
the smallness of lel/A^ and b/£ . Therefore we expect 
that at most 



r - o ' TnEp 



(4.23) 



for (p x — ip y ) superconductors with the test potential 

Km. 



The result (4.23) is unchanged even away from the 
Born limit within the t-matrix approximation. Crucial 
is the structure of M_ in £; it is common for Born and 
t-matrix approximations. 



A comparison between ( 3.29 ) and ( 4.25 ) clearly shows 
that the impurity effect in vortex of superconductors 
with time-reversal-symmetry-breaking depends severely 
on the chirality for Eg <C E v <C A^ even under the same 
trial amplitude Ao(r) and vorticity of the pair-potential. 

In the vortex core with +1 vorticity of (p x — ip y ) su- 
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perconductors, the result in ref. O seems to be valid in 
wider energy region while the random-matrix approach 
becomes less applicable, compared to the (p x + ip y ) case 
with +1 vorticity. 



§5. Implication for Flux Flow Ohmic Conduc- 
tivity 

Now we consider the implications of our results to the 
flux flow conductivity. The temperature dependence of 
t he flux flow conductivity in moderately clean regime 
(1.2) is governed mai nly b y T-dependences of T(E P = T) 
and Ap(r) (or £i in ( |2.1fj| )). The ratio £i/£o is expected 
to beEl* ~ T/Aqo (Kramer-Pesch effect) if the zero en- 
ergy Andreev bound states exist on the quasiclassical 
trajectory with vanishing impact parameter £3' Here we 
assume that the Kramer-Pesch effect occurs. Under this 



d@© 



assumption, the expression (1.3) reduces to 

Aoo \ H C 2T n 
~T~ j BT(E P = T) 



at/a n ~ In — - 



In 



Aqc 

T 



BT(E p = T) ■ 



(5.1) 



For the Born limit of s-wave case, the substitution, of 
Q3.29D into (|t|) reproduces the known result (|TTT|) M Ll ) 
For p-wave with d = z (p x + ip y ), we obtain from (5.1) 
and 04.160 



0(1) 

(E S /T) 2/3 



for E s « T « A^ 
for Tl/Ej < T < Eg. { ' ' 

Thus we expect novel enhancement <7f at temperatures 
lower than Eg for the present case. 

For p-wave with d = z (p x — ip y ), on the other hand, 
the impurity scattering rate T is much smaller than T n 
in energy region E < A m . The flux flow conductiv- 
ity <7f is, therefore, expected to be much larger than 
<r n [In (Aoo/T)] (H C 2/B) at temperatures much lower 
than Aqq. 

§6. Discussion 

So far we have seen that the impurity scattering rate 
of the Andreev bound states in vortex core depends sen- 
sitively on the parity and chirality of pair-potential. We 
discuss this point in more intuitive way in the present 
section. 

The impurity scattering rate of an in-going state is, in 
general, detemined by two factors: the density of states 
of available out-going states and the matrix element. We 
considered s- and p- wave superconductors having com- 
mon pair-potential amplitude (Ao(r)) and spectrum of 
the Andreev bound states. Therefore, the phase-sensitive 
scattering rate should originate from the matrix element. 
In the course of calculations, actually, we note that dif- 
ferent structures of matrix parts (3.7) (4.4) ( 4.1S| ) of the 
green function are crucial factors giving different scat- 
tering rates. Difference among the three matrices lies in 
the off-diagonal element, which is nothing but the "Ric- 
cati amplitude " 7 or 7^ of zero energy bound states on 
the trajectory at zero impact parameter in pure super- 



conductors. It has been mentionecSE^ that the 7 and 
7T in pure superconductors correspond to the ratio of 
particle amplitude u to the hole one v of the Andreev 
equation. EP In low energy limit, the amplitude of the 
two are the same and the relative phase is fixed strictly 
by the phase of the pair-potential so that the multiple 
Andreev reflections within the core constructively inter- 
fere. The phase of the pair-potential along the trajectory 
with zero impact parameter, of course, depends on the 
parity of pairing symmetry. Let us consider the scatter- 
ing rate due to Born scatterers of a bound state with 
energy e with impact parameter b e = e/E'. The scat- 
tering rate of a bound state is given by the integration 
of the scattering rate at each point along the trajectory. 
The density of states of available out-going state is given 
by the local density of states ((e — E p + i<5) _1 ), which is 
proportional to 1/s for b e <C s <C £o- This factor alone 
gives r/r n ~ ln(£o-E'/e). Thi s expectation reproduces 
the result in s-wave case ( |3.29| ) and do not give the re- 
sults in p-wave cases. Owing to the long tail of the local 
density of states as a function of s, the main contribu- 
tion to the scattering rate comes from the large-s region, 
where the scattering angle is near 7r; backward scatter- 
ing. From the argument of the previous paragraph, it is 
understandable that the strength of the backward scat- 
tering between the Andreev bound states in vortex core 
is parity-dependent. 

Chirality dependence of the scattering rates in the 
Born limit comes from the small s region (forward scat- 
tering region). The phase of the pair-potential has con- 
tributions from vorticity and chirality. Owing to the cir- 
cular symmetry in the present study, the completely de- 
structive interference occurs in the matrix element in the 
case of A(r) = Ao^"^. 

Next we consider the non-Born case of A(r) = 
A e i( * +Q) . The factor (/)(/) - (g) 2 in the denominator 
of (2.S) behaves as (£,o/s) 3 . Therefore the main effect of 
|<$o| is the suppression of the scattering rate in the region 
b € < s <C so = |<5o| 2 / 3 £o- This fact explains why the 
scattering rate decreases appreciably and chirality de- 
pendence becomes small in the energy region lower than 
E s . 

Lastly, we make a remark on the formulation used in 
§ 3 and § 4. The impurity scattering rates of the Andreev 
bound states within s-wave vortex core was already cal- 
culated in ref. Eq in a different method. Application of 
the method used in ref. E9 to p-wave case reproduces the 
result in §4. In this sense, the two methods are equiv- 
alent. The meaning of the approximation is, however, 
clearer in the present work than in the earlier work. 

The method in ref. Eq was a generalization of that of 
Kramer and Peschjj^ which is originally developed for 
pure superconductors. The spirit of the crucial assump- 
tion in the original Kramer-Pesch paper (which is writ- 
ten in the sentense immediately below eq. (14) in ref. 0) 
is not obvious as it stands. The method developed in 
the present paper is, on the other hand, a generaliza- 
tion of that presented in Appendices G and H of ref. |32| 
to impure superconductors. While the Green function 
becomes singular near the bound state pole, Riccati am- 
plitudes 7, 7^ are regular as functions of energy e and 
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impact parameter b. That is why we can expand those 
amplitudes safely with e and b and the spirit and valid- 
ity of approximation are clear. This "good" behaviour 
of Riccati amplitude comes from the behaviour of wave 
functions; we recall that the Riccati amplitudes are re- 
lated directly to the wave functions of Andreev equation 
in pure case. Of course, wave function cannot be defined 
any more after the impurity-average. The Reccati ampli- 
tudes are, however, well-defined quantities analogous to 
wavefunctions even in the presence of energy dissipation. 
Riccati formalism of quasiclassical theory of supercon- 
ductivity turns out to be powerful in analytical as well 
as numerical study. 

§7. Conclusion 

In summary, we presented the calculations of the im- 
purity scattering rate of Andreev bound states within 2D 
pancake vortex of s-wave and chiral p-wave moderately 
clean superconductors and further considered flux flow 
conductivity. 

We found that energy dependence of impurity scatter- 
ing rate is governed by the relative phase between the 
particle and hole channels in the wavefunction of an An- 
dreev bound state. This relative phase is strictly fixed 
by the parity and chirality of pair-potential and thus 
the scattering rate and resulting flux flow conductivity 
become parity- and chirality-sensitive in the moderately 
clean unconventional superconductors. The calculation 
of the impurity scattering rate under self-consistently de- 
termined pair-potential and self-energy is left as a future 
problem. 
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Appendix: Derivation of ( 4 . 1 6j ) 

In this appendix, we calculate the impurity scattering 
rate of Andreev bound states in vortex core (with +1 vor- 
ticity) of p-wave superconductors with d = z (p x + vp y ) 
within the t-matrix approximation. We consider the case 
where | <5"o | ^ 1; we hence replace cos<5o by unity and 
sin(5o by <5o in the following part. 

Following the same proceduce as that used in §4.2, we 
obtain 



g ~ nve 



- u ^M+[2C(e-E p + S)]' 



with 



s = C 



ds£(s)e" 



-u(s) 



(AT) 



(A-2) 



The imaginary part of S yields the impurity scattering 



rate T. Further the expression (A-2) reduces to 



r„/ (E s /e) . 



Here /(A) is defined, for complex A satisfying RcA > 
and ImA < 0, as 



/(A) 



-f-ooA 



Az- 



C(z) 



1 



vA c(*) (C(*0 + 1) - A a 

where C(z) is given by 



(A-3) 



i((z) = \/z 2 -l 

12 i|l/2 



If " exp [i {Arg(z - 1) + Arg(z + 1)} /2]A4) 

Principal values of arguments are, respectively, taken as 

< Arg(z - 1) < 2tt, -it < Arg(z + 1) < tt. (A-5) 

Corresponding branch cuts are given by Imz = and 
|Rez| > 1. The path of integral in (A-3) is taken within 
the Riemann surface specified by ( |A-5[ ). From now on, 
we impose the condition that — 7r/4 < argA < 0, which 
is equivalent to < arge < 7r/4 and evaluate the integral 
in the tw o lim its |A| 3> 1 and |A| <§; 1. In both limits, the 
integral (A-3) can be written as the sum of the residue of 
the pole of the integrand and the integral along a branch 
cut; 



/(A) =ipole(A)+Icont(A). 

The contribution from the pole zq which yields 



C(*o) (C(*o) + 1) = A 2 



(A-6) 



turns out to be 



^polc(A) — 



27riC(z ) 



z (3C(«o) + 1) 
The other contribution 7 con t is given by 

y 2 (y 2 + 



t (A) 



-2i 



dy- 



1 + A 2 



V / ^TT{(y 3 -y) 2 + (2y2 + A2) 2 } 

(A-7) 
First we consider the asymtotics of / po io(A). In the 
limit |A| < 1, it follows from ( |A-4| ) and QA-6| ) that 
C(^o) - * A 2 and zq — * 1. We immediately see that 



7 pole (A) -» 2mA 2 for |A| 



0. 



(A- 



When |A| ^S> 1, on the other hand, we obtain that 
C( z o) —* A 2 / 3 and zq — > iA 2 / 3 . Here the branch cut in 
A plane is taken along the real axis with ReA < 0. These 
results lead to 



/poie(A) -► 2^A- 2 / 3 /3, for |A| 



(A-9) 



Next consider the behavior of 7 con t(A) in the two limits. 
Around A = 0, icont(A) is a regular function of A and is 
expanded as 

44iA 2 
7 cnnt (A) = -2i+^^+0(A 4 ). (A-10) 
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When 



A| 3> 1, on the other hand, the main contribution 
in (A-7) comes from the region j/>1. We can, therefore, 
ignore unity in the integrand and ( |A-7| ) reduces to 

dr? (tj + A 2 ) 



dr^- 



(ATI) 



r] (77 - 1) + 4t7 2 + 4A 2 ?? + A 4 
The integrand in ([ATI) approaches A 2 / (r; 3 + A 4 ) for 
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77 < |A| 2 and 1/(77 + 4A 2 ) for |A| 2 < 77. The former 



contribution to (A- 11) turns out to be 



A 2 d?7 



rf 



A 2 d?/ 



vr + 



A 4 



27TJA- 2 / 3 
3V3 



(A-12) 

where the branch cut in A plane is taken along the real 
axis with ReA < 0. On the other land, the contribu- 
tion from the region |A| 2 -C 77 can be easily shown to be 
C(l/A). Therefore ( A-12 ) gives the asymptotics of 7 cont 
in the limit of |A| — > 00. 

Now we find from ([a|) and ( |A-10|) that S in ( |I^ ) 
becomes 

- 2if n (1 + 0((e/E 5 ) 2 )) , for |e| » Eg (A-13) 

and from (O) and (|A-12|) that 



2tt/3 (l - i/V3) f n (e/E s ) 2/3 , for |e| « E 5 . (A-14) 
Expression (|AT3|) immediately leads to T ~ 2r n ; re 



suits in (4.16) for that energy region. In (AT4), on the 
other hand, a real part exists, which imposes additional 
restriction on the validity of our results. We have as- 
sumed that the dispersion of the bound stat e en ergy is 
proportional to the impact parameter b (see (3^) in §3). 
This assumption holds only if the condition 

e»f n (e/%) 2/3 (A-15) 

is satisfied. We note that e ~ E p in the vici nity of the 
bound state pole . Consequently, we obtain ( 4.16 ) from 
(|A-14|) and ( |A-15| ). 
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